Linear "ship waves" generated in stationary flow of a Bose-Einstein condensate past 

an obstacle 
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Using stationary solutions of the linearized two-dimensional Gross-Pitaevskii equation, we describe 
the "ship wave" pattern occurring in the supersonic flow of a Bose-Einstein condensate past an 
obstacle. It is shown that these "ship waves" are generated outside the Mach cone. The developed 
analytical theory is confirmed by numerical simulations of the flow past body problem in the frame 
of the full non-stationary Gross-Pitaevskii equation. 
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INTRODUCTION 

Experimental creation of Bose-Einstein condensate 
(BEC) has led to emergence of a new field of nonlinear 
wave dynamics owing to a remarkable richness of non- 
linear wave patterns supported by this medium. First, 
vortices and bright and dark solitons were observed and 
their dynamics was studied theoretically in framework 
of the mean-field approach (see, e.g. [l[ and references 
therein). Then, dispersive shocks generated by a large 
and sharp disturbance of BEC were found in experiment 
0, 01 and explained theoretically [H, Q in framework of 
the Whitham theory of modulations of nonlinear waves 
(see also numerical experiment in JBj]). At last, the sta- 
tionary waves generated by supersonic flow of BEC past 
obstacles have been recently observed [6j. They were 
studied in 0, H, Q where the main focus was on the non- 
linear component representing a train of solitons (in the 
simplest case of a single soliton Q) or, more precisely, 
having a form of a modulated nonlinear periodic wave. 
The theory developed in 01i shows that there exist 
stationary spatial solutions of the Gross-Pitaevskii (GP) 
equation which describe nonlinear waves supported by a 
supersonic BEC flow with Mach number 

u , „ 

M = - > 1, 1 

c s 

where u is velocity of the oncoming flow at x — > — oo 
and c s is the sound speed of the long-wavelength linear 
waves. The density n of the condensate, as well as the 
components of the velocity field, depend on the variable 

w = x — ay (2) 

alone, where a is the slope of the phase lines with respect 
to the y axis and it is supposed that the velocity of the 
oncoming flow is directed along the x axis. Then, the 
Mach cone for sound waves with infinitely large wave- 
length corresponds to the slope 

a M = VM 2 - 1, (3) 



and it was shown in [9J] that the spatial (oblique) solitons 
have a > a^.j, that is they are located inside the Mach 
cone. In particular, it was shown that the shallow solitons 
are formed close to the Mach cone, a — au <C clm, and 
they are asymptotically described by the Korteweg-de 
Vries (KdV) equation; and deep solitons have a^S> 1 (i.e. 
they are formed at small angles with respect to direction 
of the oncoming flow) and are asymptotically described 
by the nonlinear Schrodinger (NLS) equation. On the 
contrary, the linear waves are generated outside the Mach 
cone and they have a < clm- In fact, these linear waves 
had been observed in numerical simulations some years 
ago [l0( but the theory of their generation has not been 
developed so far. The aim of this paper is to develop 
such a theory. 

LINEAR WAVES GENERATED IN A BEC FLOW 
PAST AN OBSTACLE 

Our analysis is based on the use of the mean-filed de- 
scription of BEC dynamics in the framework of the GP 
equation 

at Ira 

where ip( r ) is t ne order parameter ("condensate wave 
function"), U(r) is the potential which confines atoms 
of a Bose gas in a trap and/or describes interaction of 
the BEC with the obstacle, and g is an effective coupling 
constant arising due to inter-atomic collisions with the 
s-wave scattering length a s , 

g = 47r?i 2 a s /m, (5) 

m being the atomic mass. Here we consider the Bose- 
Einstein condensate with repulsive interaction between 
particles for which g > 0. 

As suggested by the experimental set-up jf^, we con- 
sider a two-dimensional flow of a condensate, so that the 
condensate wave function ip depends on only two spatial 



2 



coordinates, r = (x, y). To simplify the theory, we as- 
sume that the characteristic size of the obstacle is much 
less than its distance from the center of the trap, so that 
the oncoming flow can be considered as uniform with 
constant density no of atoms and constant velocity Uo 
directed parallel to x axis (see also estimates in Q). It 
is convenient to transform Eq. ([5]) to a "hydrodynamic" 
form by means of the substitution 



tp(r, t) = y/ n(r, t) exp ( ~ J u(r',t)dr' 



(6) 



where n(r,i) is density of atoms in BEC and u(r, t) de- 
notes its velocity field, and to introduce dimensionlcss 
variables 



y = y/V2£, t = {c s /2y/2^% 



n/no, 



u/c s 



(7) 



where £ = h/ y/2mnog is the BEC healing length and nu- 
merical constants are introduced for future convenience. 
As a result of this transformation we obtain the system 
(we omit tildes for convenience of the notation) 



|u t 



(uV)u + Vn + V 



(Vn) 2 
8n 2 



V(nu 
An" 

An 



(8) 



where V = (d x ,d y ). Since we shall consider waves far 
enough from the obstacle, the potential is omitted in ([8]). 

We are interested in linear waves propagating on the 
background flow with n = 1, u = M, v = 0. Hence, we 
introduce 



n = 1 + ni, u = M + m, 



v = Vi , 



(9) 



and linearize the system (j5J) with respect to small devia- 
tions n±,ui,v±. As a result we obtain the system 

\n lit + u\ <x + Mm tX + v x ,y = 0, 
|ui,t + Mui <x + ni. x - \(ni tXXX + n liXyy ) = 0, (10) 
|ui,t -I- Mvi, x + m >y - \{ni. xxy + ni my ) = 0, 

which describes propagation of linear waves in BEC with 
a uniform flow. We obtain the applicability condition of 
these equations, if we notice that in the linear wave u\ ~ 
Mn\ and the nonlinear terms of the order of magnitude 
UiVui ~ V(Miti) 2 can be neglected as long as they are 
much less than the linear ones ~ Vni. Thus, we get the 
criterion 



m < l/M 2 



(ii) 



Hence, if M is large enough, the linear theory is appli- 
cable to description of waves outside the Mach cone far 
enough from the obstacle, where the density amplitude 
ni of the wave satisfies the condition pip . For harmonic 



waves tii, ui,Vi oc exp[i(k x x+k y y) — iujt] the system (fpQ|) 
yields at once the dispersion relation 



-=Mk x ±k^l + -, 



(12) 



where k = y fc 2 + fc 2 . 

Now we consider the stationary wave patterns far 
enough from the obstacle where the condition pip is sup- 
posed to be fulfilled. In fact, this problem is analogous 
to the Kelvin theory of ship waves generated by a ship 
moving in a deep water, but with a different dispersion 
law (jX2j) . Hence, we shall follow Kelvin's method in its 
form presented in 



111, l!2| 



First, we notice that in a stationary wave pattern lj 
i) and, hence, the components of the wavevector k = 
(k x ,k v ) are the functions of the space coordinates (x,y) 
connected with each other by the relationship 



G(k x ,k y ) = Mk x + k^l + — =0, (13) 

where we have taken into account that for chosen geom- 
etry of the BEC flow the wave must propagate upwind 
i.e. k x < 0. 

Next, the "ship wave" pattern corresponds to a mod- 
ulated two-dimensional wave where the wavevector k is 



a gradient of the phase [11 1 



= / k-dr. 

Jo 



(14) 



Hence, the components (k x ,k y ) satisfy the condition 
dk x dk 



dy 



dx 



= 0, 



(15) 



which, with an account of (|13p . yields the equation for 



dk y 
dx 



f'(ky 



dk y 
dy 



0. 



(16) 



where f'(k y ) is defined by the derivative of an implicit 
function (fl3l): 



/' = - 



dG/dky 
dG/dk x ' 



(17) 



It follows from Eq. (fT6)) that k x and k y are constant along 
the characteristics defined as solutions of the equation 



(18) 



At last, since at large distances from the obstacle, the 
latter can be considered as a point-like source of waves, 



3 




FIG. 1: Sketch of the wave front with wavevector k normal 
to it at the point A and coordinates x an d V defining the 
parameters of the wave pattern at this point. 



k 2 

G{k,rj) = —Mkcosrj + kd 1 + — = 0. (21) 

Then elementary calculation gives for (|19|) the expression 

(2M 2 - 1 - tan 2 n) tan n 

tanY = — K — r 5 ^—z — '— (22) 

A (A// 2 + 1) tan 2 7] - (M 2 - 1) V ' 

and Eq. fl2T]) yields 

fc = 2^M 2 cos 2 77- 1. (23) 

Thus, we have found that for a fixed value of 77 the compo- 
nents {k x , k v ) are constant along the line \ — const with 
X defined by Eq. (|22|) and the length of the wavevector 
given by (j2"3"|) . Therefore the phase (fT4"|) can be conve- 
niently calculated by integration along the line \ = const 
with constant vector k, so that 



= (k cos [i)r. 



(24) 




FIG. 2: Wave pattern of stationary linear waves generated in 
the flow of BEC past a point-like obstacle. Dashed line cor- 
responds to the Mach cone of linear waves in the long wave- 
length limit. 



the resulting stationary wave is centered, that is the char- 
acteristics intersect at the point (x, y) = (0, 0) of the lo- 
cation of the obstacle. Hence, we obtain the solution 



V 
x 



tanx 



dG/dk y 
~dG/dk x ' 



(19) 



where x is the angle between the x axis and the radius- 
vector r of the point A with wavevector {k x ,k y ); see 
Fig. 1 where other convenient parameters are also de- 
fined, namely, the angle rj between the wavevector k 
and a negative direction of the x axis, and the angle 
/i = 7T — x ~ V between vectors k and r. 
Thus, we have 



-fccosry, k y = ksinr), 



(20) 



It means that the lines of constant phase (e.g. the wave 
crests) 9 are determined in parametrical form by Eq. (|22|) 
and 



(25) 



k cos \i ' 



where k is given by Eq. (f2"3"|) and /1 can be calculated 
from tan/x = — tan(% + 77) which gives, after elementary 
algebra, the expression 



tan fi 



2M 2 

-— sin 277. 



This expression permits one to express Eq. (J25J) as 



4:9 

— ^/M 2 {M 2 - 2) cos 2 77 + I 



and Eq. {25]) as 



tanx 



(1 + fc 2 /2)taii77 
Af 2 -(l + fc 2 /2)' 



(26) 



(27) 



(28) 



At last, the curves with constant phase 9 are given in 
Cartesian coordinates by the formulas 



V 



AO 

rcosx — 73- cos 77(1 — M 2 cos27?), 

r sinx = — $int](2M cos 2 7] — 1). 
k 6 



(29) 



Thus, we have found the expressions describing the linear 
wave pattern in a parametric form where the parameter 
?7 changes in the interval 



1 1 

arccos < ?/ < arccos 



M 



M 



(30) 



For rj — we have x < and y = 0, that is small values 
of 77 correspond to the wave before the obstacle. The 
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boundary values r\ 
lines 

x 
V 



±arccos(l/M) correspond to the 



1 = ±a M , 



(31) 



that is the curves of constant phase become asymptoti- 
cally the straight lines parallel to the Mach cone. The 
general pattern is shown in Fig. 2. 



NUMERICAL SIMULATIONS 

We have compared the above approximate analytical 
theory with the exact numerical solution of the GP equa- 
tion, which in non-dimensional units (JT]) takes the stan- 
dard form 

#t = + ^ yv ) + U{x,y)if) + HV, (32) 

which corresponds for U — to the system © with 



V^exp yi J u(r',i)dr' 



In our simulations the obstacle was modeled by an impen- 
etrable disk with radius r = 1. Such an obstacle intro- 
duces large enough perturbation into the flow to generate 
an oblique solitons pair behind it (see Q). We assume 
that at the initial moment t — there is no disturbance 
in the condensate, so that it is described by the plane 
wave function 

H x >y)\t=o = exp(iMx) 

corresponding to a uniform condensate flow. For large 
enough evolution time, the wave pattern around the ob- 
stacle tends to a stationary structure. An example of 
such a structure for M — 2 is shown in Fig. 2. A dashed 
line corresponds to the analytic theory developed in the 
preceding Section; as we see it agrees very well with the 
numerical results for M = 2. 

The condition pip indicates that the nonlinear effects 
grow up with increase of M. To demonstrate this explic- 
itly, we have compared the wavelength A at y = calcu- 
lated using the developed linear analytic theory with the 
same parameter obtained from our full numerical simu- 
lations. According to linear theory, the wavelength at 
y = (i.e. rj = 0) is constant and equal to 

(33) 



A 



1 



where we have used Eq. (f2"3")) . In Fig. 4 we compare this 
dependence of the wavelength A on the Mach number 
M with the results of numerical simulations at the point 
with ni w 0.1. As we see, Eq. (|33[) is very accurate 
for values of M satisfying the condition (fTT|) and discrep- 
ancy between analytical and numerical results slightly in- 
creases with increase of M. In general, this plot confirms 
validity of a linear theory in the region of its applicability. 



>. 




FIG. 3: Numerically calculated wave pattern of stationary 
linear waves generated in the flow of BEC past an obstacle. 
The Mach number is equal to M = 2 and the radius of an 
impenetrable obstacle to r = 1. Dashed line corresponds to 
the linear analytical theory for the line of constant phase. It 
corresponds to Eqs. 1|29[) but the curve is shifted to 1 unit of 
length to the left from the center of the obstacle for better 
fitting to numerics. It should be noted that this distance is 
negligibly small in our theory corresponding to a point-like 
obstacle. A pair of oblique dark lines behind the obstacle 
correspond to spatial solitons studied in 



aralytical 
• numerical 
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FIG. 4: Dependence of the wavelength of the wave generated 
before the obstacle on the Mach number calculated numeri- 
cally (full time-dependent GP equation) — dots, and analyti- 
cally (linear theory) — dashed line. 



CONCLUSION 

We have developed here the theory of linear waves gen- 
erated by the flow of BEC past an obstacle. The linear 
approximation is correct for small enough amplitudes of 
the perturbation. This condition is satisfied in the case 
of small disturbance introduced by the obstacle and not 
too high values of the Mach number. Our numerical sim- 
ulations confirm the analytical theory in the region of its 
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applicability. 
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